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^ . Abstract - The aim of this paper is to study the spatial behaviour of the so- 

\0 '. lutions to the boundary-final value problems associated with the linear theory 

^ of elastic materials with voids. More precisely the present study is devoted to 

porous materials with a memory effect for the intrinsic equilibrated body forces. 

An appropriate time-weighted volume measure is associated with the backward 
. in time thermoelastic processes. 

Q^c Then, a first-order partial differential inequality in terms of such measure is 

established and further is shown how it implies the spatial exponential decay of 

^— > ■ 

the thermoelastic process in question. 

1 Introduction 

X 

The boundary-final value problems associated with the linear thermoelasticity have been 
studied by Ames and Payne [1] in connection with the continuous dependence of the ther- 
moelastic processes backward in time with respect to the final data. It is well known that 
this is an improperly posed problem. A further study on this subject was recently developed 
by Ciarletta [2]. 

More recently, Chirita and Ciarletta have developed an exhaustive description for the 
spatial behaviour of solutions in linear thermoelasticity forward in time [3, 4, 5]. In this 
connection some appropriate time-weighted surface power functions are introduced. On this 
basis some spatial decay estimates of Saint-Venant type are established for bounded bodies; 
while for unbounded bodies some alternatives of Phragmen-Lindelof type are obtained. In 
the present paper, we aim to extend the above context to porous materials for backward in 
time processes. We shall refer to the well-known theory by Nunziato and Cowin, in which 
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and under auspices of G. N. F. M. of the Italian Research Council (C. N. R.). 
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the presence of small pores (or voids) in the conventional continuum model is introduced 
by assigning an additional degree of freedom, namely, the fraction of elementary volume 
that is possibly found void of matter [6, 7, 8]. Following a previous paper by Goodman and 
Cowin on granular (flowing) materials, the bulk mass density of the material is represented 
as a product of two fields, the void volume fraction and the mass density of the matrix 
material [6]. Starting from the fundamental work of Iesan in [9], we consider the porous 
materials with the memory effects for the intrinsic equilibrated body forces and we associate 
with such a model the boundary-final value problem as in [8] and [10]. Then, we study the 
spatial behaviour of the thermoelastic processes backward in time by using an appropriate 
time-weighted measure. 

In Section 2, we present the basic equations of the linear dynamic theory of porous body 
developed in [8] and [10]. Some constitutive assumptions and other useful results are also 
presented. In Section 3 a time-weighted volume measure is defined and a first-order partial 
differential inequality is established in terms of such measure. On this basis a spatial decay 
estimate of Saint-Venant type is proved; our decay results look quite similar to those estab- 
lished by Lin and Payne for the backward heat equation defined on a semi-infinite cylinder 
[11]. 



2 Preliminaries 

We shall denote by B the smooth domain of the physical space (= IR 3 ) occupied by an 
anisotropic, homogeneous and porous body in a fixed, natural reference configuration. Iden- 
tified IR 3 with the associated vector space, an orthonormal frame of reference is introduced. 
The vectors and tensors will have components denoted by Latin subscripts (ranging over 
{1,2,3}, unless otherwise specified). Summation over repeated subscripts and other typical 
conventions for differential operations are implied, such as a superposed dot or a comma fol- 
lowed by a subscript to denote partial derivative with respect to time or the corresponding 
coordinate. Occasionally, we shall use bold-face character and typical notations for vectors 
and operations upon them. 

We study the boundary-final value problems associated with the linear thermoelasticity. 
Following [8, 10], the local balance equations become 

Sjij + pfi = piii, balance of momentum, 

hi t i + g + p£ = px<P, balance of equilibrated stress, (1) 

p9 r] = q iti + pr, in B x (— oo, 0) energy equation. 

In these equations, u is the displacement vector fields; is the change in volume fraction 
starting from the reference configuration; 6 is the temperature variation from the uniform 
reference temperature #o(> 0). Moreover, S and f are the stress tensor and body force, 
respectively; h, g and t are the equilibrated stress vector, intrinsic and extrinsic equilibrated 
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body force, respectively; rj, q and r are the specific entropy, the heat flux vector and the ex- 
trinsic heat supply, respectively. Finally, p and X are the bulk mass density and equilibrated 
inertia in the reference state, respectively. We assume that f, £ and r are the continuous 
functions onBx (— oo,0], with B closure of B. 

By denoting U = {u, ip, 6}, it follows the strain fields are 

e ij = \( u i,j + u i,i)i 7» = <P,i, K i = 0,» in B x (-00, 0]. (2) 

The surface tractions s, the surface equilibrated stress h and the boundary heat flux q are 

Si = SjiUj, h = hjUj, q = qjfij, (3) 

where n is the outward unit normal vector to the boundary surface. 

We consider the porous bodies with the memory effect for the intrinsic equilibrated body 
forces. Our attention is on the materials, which are initially free from stress and the intrinsic 
equilibrated body force, the entropy and the heat flux rate are equal zero. The constitutive 
equations are given by 

hi = D rsi e rs + Aij-fj + biip - a { d, 



g = -np + G, 

G = -B^j - bdi -£(p + m9, 



(4) 



pi] = M H e k i + rrup + + aO, 

qi = KijKj, in B x (— 00, 0]. 

The material coefficients satisfy the relations 

Gijrs G rs ij Gji rS ] -Dijr -^jir^ ^-ij -^-jii ^jii -^ij I ji • (5) 

and 

r > 0. (6) 

Throughout this article, we assume the bulk mass density p, the equilibrated inertia x an d 
the constant heat a are strictly positive. 

The constant conductivity tensor K is a symmetric positive definite tensor; thus, there 
exist the positive constants k m and ku such that 

k m KiKi < KijKiKj < k^jKiKi Vk.. (7) 

The constants k m and ku are the minimum and maximum conductivity moduli for K. 
It follows from Schwarz's inequality and the last inequality that 

ftft = KijKjqi < (KijKiKj) 1 ^ (K rs q r q s ) 1/2 < {Ki j KiK j ) 1,2 {k M q s q s ) 1,<1 , (8) 
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so that 

q%qi < k M KijKiKj. (9) 

For what follows, it is useful to introduce the vector space £ of all vector fields of the form 

E = {Ei^xm,^}, with Eij = Eji, and Xi = Vx- (10) 

Moreover, for each E e £ we define the vector field S(E) as 

S(E) = |^, Xl 

where 

Cij rs E rs -\- Dijgiig -\- JBijijj^ 
hi = D rsi E rs + AijiTj + hip, (11) 
G = -BijEij - bi-Ki - £ip, 

and the coefficients obey the symmetry relations (5). The vector field S(E) belongs to £, 
too. 

Now, for any E, E G £, we consider the following bilinear form 

2.F(E, E) = C ljrs E tj E rs + ^ + AijTTiTtj + B^E^i, + 

(12) 

+D ijs (Eij7r s + EijTTs) + bi(^7Ti + 1p7Ti), 

where E = {Eij,Xi^i,i>}- From (11) we have 

2JF(E, E) = [SjiEij + hm - Gifj} and 2JF(E, E) = [S Jt E tj + km - Gip]. (13) 
By the constitutive equations, we prove 

JF(E, E) = JF(E, E), VE,Ee£. (14) 

The Cauchy-Schwarz's inequality implies 

<{W(E)] 1/2 {W(E)] 1/2 , VE,Ee£, (15) 

where W is the quadratic form associated to T 

2#(E) = 2.F(E, E) = C tjrs E tj E rs + ^ 2 + A^ttj + 2B lj ^E lj + 

(16) 

+2D ijs E ij 7r s + Ibi^Tii. 
We assume that V4^(E) is a positive definite quadratic form and, consequently, we have 

/i m (E ij E ij + X nun + V> 2 ) < 2W(E) < /i M (E ij E t] + X mm + . (17) 
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By setting E = S(E) in (17) and (13), we prove 

2W(S(E)) < hm (sjiSji + Ui-K + G 2 ^ , (18) 

and 

+ + G 2 = ZF(E, S(E)) < 2[^(E)] 1 / 2 [#(S(E))] 1 / 2 < 

X , .1/2 (19) 

< 2[#(E)] 1 /Vm \SjiSji + -ink + G 2 j ; 

so that this inequality implies 

SjiSji + -kk + G 2 < 2fi M W(E). (20) 

If E = {eij, Xilh V 9 } i then we introduce the following notations 

2W* = 2W{E) = C l3rs e i3 e ra + & 2 + A lj7i7j + 2B ij ipe ij + 2D ijs e ijls + 2b^ lu 

G = -Bijdj - k% - 
By eqs. (11), (13) and (21), we can see that 

2W* = [Sjidj + kli - G<p] , W* = [Sjiiij + kii - Gcp] . (22) 

For any positive number e and each second-order tensors, L and F, we have the inequality 

(Lij + F tJ )(L tJ + Fij) < (1 + + (1 + -jF tJ F tJ . (23) 
By aid of eqs. (4), (21) and (23), we obtain 

S tj = Sij - MijO, K = ~hi- OiO, G = G + m6, (24) 

and 



1 _ l _ 

SijSij + -hihi = (S^ - MijO^Sij - M^e) + -(h* - a^)(^ - a^) < 

(1 + e)SijSij + (1 + ^M^MijO 2 + + e)kk + (1 + -)<H<H^\ < 

6 X. £ 



(25) 



<{l + e)2 m W* + {l + \)M 2 6 2 J Ve>0, 

in which 

M 2 = max(M ii M ij + -a^). (26) 

B X 



We consider the boundary- final value problem V defined by the equations of motion (1), 
the geometrical equations (2) and the constitutive equations (4) and the following final- 
boundary conditions 



Ui(x, 0) = «?(x), ^(x, 0) = m?(x), </?(x, 0) = y9°(x) 
p(x,0) = <p°(x), 0(x,O) = 0°(x), xGB, 



(27) 



and 



ui = u* on Si x (— oo, 0], Si — s* on E 2 x (— oo, 0], 

<p — ip* on E 3 x (— oo, 0], h = h* on E 4 x (-oo, 0], (28) 

= 0* on E 5 x (-oo,0], g = q* on E 6 x (-oo, 0], 
where (i — 1, ... ,6) are the subsets of dB such that 

Ei U E 2 = E 3 U E 4 = E 5 U E 6 = dB, E : n E 2 = E 3 n E 4 = E 5 n E 6 = 0. 

The terms on the right-hand of eqs. (27) and (28) are prescribed continuous functions. 

For further convenience we use an appropriate change of time variable and notations to 
transform the boundary-final value problem V in the boundary-initial value problem V* 
defined by the following equations 

hi,i + 9 + p£ = PXV, (29) 
-p^oV = qi,i + pr, in B x (0, +oo), 

e%j = -(u id + u jti ), 7i = (p ti , m = 0,i, infix [0, +oo), (30) 

and 

Sij Cij rs c rs ~\~ Dij s —f s -\- Bij(p AdijQ, 
hi = D rsi e rs + Aij-fj + bi(p - a { 6, 
g = rip + G, 

G = -BijBij - - £ip + m9, 
prj = M H e k i + a,fii + m<p + a0, 

qi = KijKj, in B x [0, +oo), 

with the initial conditions 

Wl (x,0) = «?(x), ^(x,0) = u?(x), </?(x,0) = /(x), 

<p(x,0) = <p°(x), 6>(x,0) = 0°(x), xGfi, 



(31) 



(32) 
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and the boundary conditions 

Ui = u* on Ei x [0, +00), Si = s* on E 2 x [0, +00), 

ip = (f* on E 3 x [0,+oo), h=h* on E 4 x [0, +00), (33) 

9 = 0* on E 5 x [0,+oo), q = q* on E 6 x [0, +00). 

We define as solution of the boundary-initial value problem V* a process ir = {u, 
e, S, if, 7, h, g, 9, k, q, t]} that satisfies eqs. (29)-(33) and 

i. Ui , ip G C 2 ' 2 (B x [0, +00)), 9 e C X '\B x [0, +00)); 

ii. tij = eji, 7,, Ki G C 1 ' 1 ^ xjO, +00)); 

iii. Stj = Sji, hi,q h g G C^{B x [0,+oo)); V G C°>\B x [0,+oo)). 

3 Saint- Venant's Principle 

Let us consider a given time T G (0, +00) and a given ( external ) data V = {f , £, r; u°, 
</?°, v? , 0°; u*, s*, ip*, h*, 9*, q*} in the problem V* . We denote by D T the support of 
the initial and boundary data, the body force, the extrinsic equilibrated body force and the 
heat supply on the time interval [0, T], i. e. the set of all xG73 such that: 

i. if xG B, then 

u?(x) ^ or u?(x) ^ 0, or p°(x) ^ 0, or y?°(x) ^ 0, or 0°(x) ^ 0, 

(34) 

or /j(x,s)^0, or £(x, s) 7^ 0, or r(x, s) 7^ for some sG [0,7]; 

ii. if x G dB, then 

u*(x, s) 7^ for some (x, s) G Ei x [0, T], 
or s*(x, s) 7^ for some (x, s) G E 2 x [0, T], 
or <^*(x, s) 7^ for some (x, s) G E 3 x [0, T], 

(35) 

or h*(x, s) 7^ for some (x, s) G E 4 x [0, T], 
or 0*(x, s) 7^ for some (x, s) G E 5 x [0, T], 
or g*(x, s) 7^ for some (x, s) G E 6 x [0, T]. 

We assume that D T is a bounded set. We consider a non-empty bounded regular region D T 
such that D T C 7}£ C 75. We note that 

i. if 7^ 79^, then we choose D T to be the smallest bounded regular region in B that 
includes Dt ; in particular, we set D T = Dt if Dt is also a regular region; 

ii. if Dt = 0, then D T may be chosen in an arbitrary way. 
Now, we mean the set D r , by 

D r = {xG75 : 7^nE(x,r) 7^ 0}, r > 0, (36) 
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where S(x, r) is the closed ball with radius r and center at x. Clearly, D T C = D Q C D r 
(r > 0). Further, we set B r = B\ D r ; we have B r . 2 D B ri , and B(ri,r 2 ) = B r . 2 \ B ri for 
r\ > T2- Let L be the diameter of Bq. The surface S r is the subsurface of dB r contained 
inside B and whose outward unit normal vector is oriented to the exterior of D r . 
In what follows we need the below lemma 

Lemma 1. Let 7r = {u, e, S, tp, 7, h,g, 9, k, q, t]} be a the solution of V* . Then, for every 
regular region P C B with regular boundary dP and for each t G [0, T], we obtain 



) + 



/o /p e " S { 2 [M(S) ^ (S) + PX ^ (S) + a6 ^ s) + 2W "* (S)] + rip2{ * s) ~ 

1 /■ 1 

+— KijKiWK^sfidvds = / e A *-[p^(i)^(i) + px<p 2 (i) + a9 2 (t)+ 

Uq J p Z 

dads+ 



+2W*(t)}dv 



JdP 



JP 



Si(s)ui(s) + h(s)ip(s) - —9(8) 
pr 



pM'Ms) + pt<p(s) - WW 



dvds+ 

- j p \\P%<% + PX¥°¥° + a9°9° + 2W*(0)]dv, 



where A is a prescribed positive parameter. 
Proof. By eqs. (29)-(31) we deduce that 



d fl 



5s 12 



: [p^( S )^(s) + pxv 2 {s) + afl 2 (s) + 2W*(s)] = rp 2 (s) + 



+^-K ij K i (s)K j (s) + pfi{s)ui{s) + p£<p(s) - ^6»(s)+ 
"0 "0 



+ 



Thus, we can see that 
d 



— | e ^I[p^( s )^( s ) + px ^{s) + a9\s) + 2W*(s)}} = 



.A 



e As ^K(sK(s) + PX0 2 (s) + a# 2 (s) + 2W*(s)] + e As ry? 2 (s) + 



1 



+e "—KijKiWKjis) + e 



As 



+e 



As 



^,( S )^( S )+/ lj ( S )^( S )-^( S ) 

"0 



pr 

pfi(s)Ui(s) + p#(s) - -z-0(s) 



+ 



(37) 



(38) 



(39) 



If we integrate this relation over P x [0,T], then we obtain the desired result with the help 
of the divergence theorem and eqs. (3). • 

For a prescribed strictly positive parameter A and for any r e [0, L], t e [0, T], we associate 
with the solution 7r the following time-weighted volume measure £(r,t) (> 0) 
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£(r,t)= f [ e Xs {^lpu l (s)u i (s)+p X <p 2 (s) + a9 2 (s) + 2W*(s)] + 

Jo J Br 2^ (4Q) 

+np 2 (s) + —KijKi(s)Kj(s)}dvds. 

The parameter A, in the above function, is characteristic for the considered measure. 
Taking into account that, for r\ > r 2 

ft r \ 

/ \ c r A i 



£(r u t) - £(r 2 ,t) = - f f e Xs {-[pu t {s)u t {s) + p X ^ 2 {s) + 

JO JB(n,r 2 ) 2 

+a6 2 (s) + 2W*{s)\ + Tip 2 (s) + ^-KijKi{s)Kj{s)}dvd8, 



Oo 

it is a simple matter to prove the following lemma. 

Lemma 2. Let 7r be a solution of initial-boundary- value problem V* and Dt be the 
bounded support of the external data T> on the time interval [0, T\. Then, the corresponding 
time-weighted volume measure satisfies the following properties 

(i) £(r,t) is a non-increasing function with respect to r, i. e. 

£(r u t) < £(r 2 ,t), with n > r 2 ,t e [0,T]. (42) 

(ii) £(r, i) is a continuous differentiable function on r e [0, L], t e [0, T] and 



<9r 



£{r,t) = -f I e Xs {-[pu,(s)u t (s)+px^(s) + a0 2 (s) + 2W*(s)] + 

JO JS r 2 



+T(p 2 (s) + —K i jK i (s)nj(s)}dads, 
u 



(43) 



-£(r,t) = jf e ^{-[p^(t)^(t) + px^ 2 (t) + a^(t) + 2W(*)]+ 

+T(p 2 (t) + —KijKi^Kj^ydv. 



(44) 



Lemma 3. Let 7r be a solution of the initial-boundary-value problem V* and D T be the 
bounded support of the external data V on the time interval [0,T]. Then, £(r,t) satisfies 
the following first-order differential inequality 

£(r,t)<-^£(r,t)+ ~S(r,t), Vr e [0, L], t e [0, T], (45) 

where 



C(A) = J^(l + e), (46) 



and 



1 M 2 AA: M 
!+£=-+ ^ + ,„ + 



3/2 ^ M y + J ^ (47) 



2 2appM 4:9oapM \ \2 2app M 4a# ^M / ap^M 
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Proof. Taking P = B r into Lemma 1, we have 

1 



.As 



£(r, () = 

-/'/ 

Jo Js, 

It follows from (44) and (48 



e^^pu^iait) + px^ 2 (t) + a6 2 (t) + 2W*(*)]*>- 



Si(t)ui(t) + h(t)<p(t) - J-0(t) 



(48) 



dadt. 



J b e Ml -[pu i (t)uM + PX^t) + ae\t) + 2W*(t)]dv< ~(r,t). (49) 

If we use Schwarz's inequality and the arithmetic-geometric mean inequality, we obtain 

1 

ISji^rijUiis) + hj(s)nj<p(s) - T 9(s)q j (s)n j \ < 



X 



Sji(s)Sji(s) 



+ 



1 r A , 2r, . 2 . £i rA/i ? (s)/i 7 (s) 

Aei 2 A A(p + ^)L 2 X 

+ ^k [ 2 a ^ 2(s)] + Br qjis)qj{s% V£l '" 2 > °- 



By eqs. (9), (25) and (50), we deduce that 



(50) 



|%(sK-Mi(s) + hj(s)nj(p(s) - —e(s)q j (s)n j \ < 



< T^K(A«i(s)ii(s) + PX^is)) + r^ 2 ( S )] + - l( ' + ( W 



+ 



Aei L 2 
£iM 2 | 1 1 
Apa e A6> £ 2 



Ap 



[AW"(a)] + 



(51) 



[\tf 2 ( S )] + ^[1^.^( S ) S .( S )], V£,£ 1; £ 2 > 0. 



If we consider £, £i, £ 2 such that 

1 ^(l + e)/^ £ X M 2 



Xp 



Xpa 



In 1 
(1 + -) + 



£ 2 & 



then 



e A6* £ 2 2a 
2aC(A) 

C(A)' 52 Xk M ' 
and e satisfies to (47), i. e. e is the positive root of the algebraic equation 



e 2 + 2e(-- 



M 2 Xk M M 2 _ Q 



2 2app M 49 ap M app M 



(52) 



(53) 



(54) 
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By eqs. (51), (53) we get for any r G [0, L] and t G [0,T] 

— / / e Xs [si(s)ui(s) + hi{s)njLp(s) — —9(s)qj(s)rij]dads < 

JO JSr t>0 /rr\ 

ft r 1 CdS ^ ' 

< I / / e Xs [si(s)ui(s) + hi{s)nj(f{s) - —9(s)q j (s)n j }dads\ < ~\^-(r,t). 
Jo J s r "o A Or 

Therefore, the relations (48), (49), (55) lead to (45). Thus, the proof is complete. • 

Theorem 1. Let n be a solution of initial-boundary- value problem V* and Dt be the 
bounded support of the external data V on the time interval [0, T]. For A (sufficiently large), 
to G [0,T], r G [0,L] such that 

L < C(A)t + < C(A)T, (56) 

we have 

S(r, t + 1^=1) < £(0, t + ^y) exp (-^r), for r G [0, L]. (57) 



Proof. We put 



l(r,t) = eM-^r)£(r,t), Vr G [0,L], i G [0,T]; (58) 

the inequality (45) takes the following form 

^(r,t)-^y^(r,t)<0, Wre[0,L],te[0,T]. (59) 

From eqs. (46) and (47) it is trivial to see £(A) ~ A 1 / 2 for A — > oo and so £(A) is an 
increasing function for sufficiently large values of A. Thus, we can choose t and r satisfying 
(56) and hence 

0<t + r ^^<T VrG[0,L]. (60) 

To — r 

By setting t = t + in (59), it follows 

C(A) 

|;[X(r,to + ^)]<0, VrG[0,L], (61) 

and therefore, 

0< X( r ,t + r -^) < J(0,t + ^y), VrG [0,L]. (62) 

The relations (58) and (62) lead to (57) and the proof is complete. • 
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